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Abstract. The heat kernel on the symmetric space of positive definite Hermitian matrices 
is used to endow the spaces of Bergman metrics of degree A: on a Riemann surface M with 
a family of probability measures depending on a choice of the background metric. Under a 
certain matrix-metric correspondence, each positive definite Hermitian matrix corresponds 
to a Kahler metric on M. The one and two point functions of the random metric are 
calculated in a variety of limits as k and t tend to infinity. In the limit when the time t 
goes to infinity the fluctuations of the random metric around the background metric are the 
same as the fluctuations of random zeros of holomorphic sections. This is due to the fact 
that the random zeros form the boundary of the space of Bergman metrics. 


1. Introduction 

In a recent series of articles M, the authors have been investigating a new approach to 
defining ‘random surfaces’. The main idea is to define integrals over the infinite dimensional 
space IC[ U0 ] of metrics of fixed area 2n in a fixed conformal class [coo] on a Riemann surface 
M as limits 

[ F(g)e~ s ^Vg := lim [ F k (g)e~ s ^V k g (1) 

JiC [ uq] k ^°° JBk 

of integrals over finite dimensional spaces B k of Bergman metrics. Given a background 
metrical uiq and a choice of a basis {sj(^)} of holomorphic sections of L k —>• M, the spaces B k 
can be identified with the non-positively curved symmetric space VN k '■= SL(N k , C)/ SU(N k ) 
of positive definite Hermitian matrices. The general question is to find sequences {d/u, k = 
e~ Sk ^V k g} of measures on B k which are independent of the choice of the basis {sj(z)}, 
which vary in a simple way under the change of the reference point oj 0 G /C[ Wo ] and have 
good asymptotic properties as k —> oo. It would be particularly interesting to construct 
a sequence {d/a k } which tends to Liouville theory measure on metrics of fixed area f M ca 0 , 
although that is not the aim of the present article. 

The sequence of measures we study in this article are the heat kernel measures 

dt&(P) Pk(t, I, P)dV(P), ( 2 ) 

where dV ( P ) is Haar measure, p k (t, P±, P 2 ) is the heat kernel of the symmetric space VN k and 
/ is the identity matrix. Under the matrix-metric identification B k ~ V]y k the identity matrix 

1 With some abuse of notation we make no distinction everywhere between the metric g and the corre¬ 
sponding Kahler form w, connected as w = ig a idz a A dz b . 
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corresponds to the background metric and the heat kernel measure ((21) is transported 
to Bk- The measure is invariant under the action of the unitary group U(Nk). Hence it is 
invariant of the choice of the basis of sections used to identify metrics and matrices. Then 
(J2J) is the probability measure on Bk induced by Brownian motion on TV*, starting at the 
identity / and continuing up to time t. The heat kernel measure is almost canonical, the 
only choices being the time t and the background metric used to make the identification 
and to start the Brownian motion. The purpose of this article is to study the behavior of 
the heat kernel measure (J2J) on Bk as k —> oo. The main geometric quantities we study are 
the area statistics 

Xu(u) = f uo (3) 

Ju 

measuring the area of an open set U C M with respect to the random area form u G Bk- We 
determine the means and variances of these random variables and their smooth analogues 
Xf(oj) = J M foj with / G C°°(M ) in various regimes, e. g. when the time t = t k is allowed 
to vary with k. The calculations are valid for any choice of background metric and the 
dependence on the background metric is simple and explicit. 

The heat kernel measure (J2]) is U (iV/J-invariant in the P variable. Such invariant mea¬ 
sures have generic form dpB k (P) ^ r n fe (e A )d/iHaar(T > ), where Ps k {& X ) is a function of the 
eigenvalues of P. It was shown in [9] that the eigenvalue density Tt 3 k (e x ) induces a function 
Tk, 2 (^ 1 , v 2 ) on M+, so that the 2 -point correlation function has the form, 

X 2 ,k{zuz 2 ) := E k 4 > p( z i)Pp(z 2 ) = Pi(zi)(j)i(z 2 ) + 7 ^/ 2 ,fc(p), (4) 

where 

\B k (z u z 2 ) | 2 

is an important invariant of the Szego kernel B k (zi,z 2 ) of the background metric, known as 
the Berezin kernel. Thus, l 2 ,k(zi, z 2 ) is the bi-potential of the variance of the area forms (or 
Kahler metrics in higher dimensions) relative to the exterior tensor product ujq Kl cno, 

Var (a;^) =E(u^ a;^) - E(a^) Kl E(a^) = E(a> 0 IE a- cu 0 E w 0 , (6) 

in the sense that 

Var (cu^) = -^(idd) z (idd) w I 2ik (z,w), . (7) 

The general formula for I 2i k{p ) for any U (V fc )-invariant measure is, 

1 f f n A + \A4 2 — B 2 

hM = A / / log(ffi 2 cos 2 0 + v\ sin 2 p) log--- p2,k( u ) sin/3 dpdu x du 2 , ( 8 ) 

z Jo z 

with 

( A = (z / 2 cos 2 p + z/f sin 2 P)p + (z / 2 sin 2 p + z/f cos 2 p) (1 - p ), 

| B — 2(z / 2 — z/ 2 )^ p( 1 — p) cos P sin p. 

The transform T Bk {e x ) —> JF 2 ,kiy) is very difficult to evaluate, and we do not know how 
to do so directly even for the heat kernel measure. The first term of (J4J) is the potential 
of uj^ I {z\)oj ( j >I {z 2 ) where a )</ )I is the background metric, and the second term h,k{p) is the 
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correction to this term, which we call the variance term. The key point is that l 2 ,k{p) is a 
function only of the variable p. This result defnes a transform 

^B fc (e A ) ->■ h,k{p) 

from eigenvalue densities to variance terms depending only on p. It would be interesting to 
know if this transform is invertible in some sense, so that one can construct U (iVfc)-invariant 
measures with prescribed pair correlation functions. In this article we calculate / 2) fc(p) when 
Fi 3 k (e x ) comes from the heat kernel measure, by a different method (also used in [9]). 


1.1. Main results. In the case of heat kernel measures we calculate the pair correlation 
function explicitly not just for Riemann surfaces, but for general projective Kahler manifolds 
(see §4.2p . The calculations give an explicit formula for the variance term of the 2-point 
function I 2k (t,p): 


2 1 e t / 2 J\ — p f°° e 21 a2 cosh A \/ coth 2 A — p + Jl — p 

- j -- / «A— = log— , =- . 

P a/27t t P J- oo a/ coth" A — p \J coth 2 A — p — a/1 — p 


(9) 


We do not integrate the result because J 2 ^(i, p) is the expected “bi-potential” and the fluctu¬ 
ations of the metric are obtained by differentiating it. An important aspect of (JUJ) is that the 
expression has no ^-dependence, except for the variable p (j5jl . which has the form e ~ kD ( z u z 2 ) 
where D(zi, z 2 ) is the diastasis (an analog of distance-squared function for Kahler manifolds) 
between the points, with respect to the background metric, see $21 Eq. (1T8|) . 

We consider several limits of this joint formula in §HJ From the geometric viewpoint, the 
most natural scaling of the time variable is t k = e^ 2 t so that the excursion distance of the 
Brownian motion in V^ k at time £*. is essentially distance t in the Mabuchi metric on /C[^ 0 ]. 

Similarly to mm, we also study the large k asymptotics both in the unsealed and scaled 
regimes. The scaling limit is common in related problems in the physics of N particles 
where one lets the number N —» oo and the volume of the surface V —» oo such that A 
tends to a limiting density (see e.g. [H] for a similar scaling in the quantum Hall effect). 
The natural length scale for metrics in the Bergman space Bk is ^=. We consider pairs 

{z\,z 2 ) G M x M to be close to the diagonal if d(zi, z 2 ) < and to be ‘off-diagonal’ if 
d{zi,z 2 ) > C where d(z\,z 2 ) is the distance between points relative to the background 
metric. In the scaling limit we consider the asymptotics of I 2 ,k(t,p) f° r pairs of points of 
the form (z,z + -^=) with \u\ < C log /c, in which case p ~ . The scaling asymptotics 

combined with the time scaling tk —> oo has a limit correlation function with a logarithmic 
singularity along the diagonal zi = z 2 (where p = 1). The variance ( idd) Zl (idd) Z2 I 2} k(t, p) 
of the Kahler metric then has a 5(zi — z 2 ) singularity along the diagonal. In fact, the scaling 
limit correlation function turns out to be identical to that for zeros of random holomorphic 
sections determined in [T7] (Lemma 3.7). 

When t G R + is finite and fixed, then / 2 ,fc(t,p) is smooth at p = 1 (see (l38jl ). and has a 
convergent expansion ( idd) Zl (idd) Z2 I 2j k(t, p) ~ \z\ — z 2 1 2 and there exist coefficients a n (t) so 
that 

OO 

h ,k(t,p) = ^ a n {t)p n . 

n =0 


(10) 
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Off the diagonal, p —> 0 and we can obtain the asymptotics by Taylor expanding the ampli¬ 
tude of the integral (]9]) • The first term ^ is singular as p —> 0. But in 1 14.31 resp. in 1 14.41 it 
is shown that the ^ ‘singularity’ cancels in the sum of the two terms. Hence the variance of 
the potential is exponentially decreasing off the diagonal. 


1.2. Comparison with [47] [18]. In 1 14.41 we first let t —> oo and then let k —>- oo. It turns 
out that in this limit, the random metrics we obtain are identical with random zero sets 
of holomorphic sections of the kth power L k of the line bundle with Chern class [wo]. As 
explained in §5] as t —> oo the mass of the heat kernel concentrates on the ideal boundary of 
the symmetric space, where the metrics correspond to the zero sets of holomorphic sections. 
On a Riemann surface, the random metrics become normalized sums of delta functions on 
random point configurations with k points. We verify that the pair correlation function of 
random Kahler metrics in B k in the limit t —> oo is given by the same formula as in mi 
for correlations between zeros of random sections. For large t (depending on k ), the random 
metric is close to such a point configuration measure. 

We now explain this similarity in more detail so that the notation and purpose of this 
article are synchronized with those of mm- The model of Gaussian random holomorphic 
sections and the results are described in more detail in 1 15.11 In those articles, the role of 
the area form (in complex dimension one) is played by the zero set measure Z s of a random 
section s G H°(M, L k ), which defines a “singular metric”. Hence it is not surprising that 
there are relations between random smooth metrics and random zero sets. 

In [47], 48], the zero current is given by Z s = ^<9<9 log |s| 2 and, analogously to ((6]) the 
variance current of zeros is defined by, 


Var(Z s /c) = E(Z s k Z s k^j — E (Zgk'j KIE(.Z s k). 


( 11 ) 


In [47], it was shown that the bi-potential I 2i k of the variancqj (fill) in the sense of ([7]) is 
given by dilogarithm 


I 2 ,k(z,w) = - 


log(l - s) 


ds . 


( 12 ) 


In §4.41 we show that the t —> oo limit of the heat kernel ensemble gives precisely the same 
bi-potential, hence the same variance. 

For zeros of random holomorphic sections over a Riemann surface, the area statistic (ED 
with respect to Z s counts the number of zeros of s in U and is denoted by Mu hr HD it is 
shown there that 

Var (Mu) = - / d zl B S2 Q k (z 1 ,z 2 ) . 

JdUxdU 

and that the number variance for zeros has the asymptotics, 


Var (Mu) = k 1/2 


i/i VoR (8U) + 0{k~^ +£ ) 


where u\ = ■ Thus, the same formula is valid in the limit f —^ oo of heat kernel random 

metrics. Similar formulae for variances of Xu and for random metrics can be derived 
from the explicit formula 0 for l 2 ,k(t, p) for any t in the heat kernel measure ensembles. 
The details are lengthy and will be presented elsewhere. 


2 In the notations of Eh h,k corresponds to 47r 2 Q/c- 
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1.3. Asymptotic central limit theorem. It is shown in [TSj that the fluctuations of the 
smooth linear statistics Xf in the case of random zeros tend to a Gaussian field with variance 
A/"(0, y/nl H 09 /H 2 ), see 1 15.21 for more details. This result holds when t — 00 for heat kernel 
random metrics. The analogous results for (j3J) do not seem to be known at present. 

It is very plausible that for general times t k the smooth linear statistics Xf with respect to 
the heat kernel measure d/j^ are also asymptotically normal, with a related variance. Such 
an asymptotic central limit theorem would be a concrete measure of how closely heat kernel 
random metrics compare to random singular metrics defined by point processes of random 
zero sets, or to those studied in (HUD Eld, where the fluctuations of linear statistics of 
eigenvalues tend to a Gaussian free field. We plan to investigate the asymptotic normality 
of fluctuations of random metrics in future work. These asymptotic normality results would 
also give a comparison of heat kernel random metrics to Liouville random metric. The 
fluctuations in the latter case are of the type of Gaussian multiplicative chaos. 

1.4. Discussion. Heat kernel random metrics are the metrics obtained by starting at the 
background metric c o ( f >I and following a Brownian motion on V^ k for time t. 

The formula (jUJ) reflects the geometry of Brownian motion of the non-positively curved 
symmetric space V ^ k , which is very different from that of Euclidean space, see f|3] for back¬ 
ground. First, due to non-isotropic nature of the Haar measure, the heat measure is concen¬ 
trated along the SU(N k )- orbit of a distinguished element h/v fc , the half-sum of the positive 
roots. Second, in the radial direction the heat kernel measure concentrates in a kind of an¬ 
nulus of radius t around the S77(iVfc)-orbit of <bv fc , see 1 13.21 This is why, as t —> 00 , the heat 
kernel measure becomes supported on the ideal boundary <9oo VN k - Its SU(N k ) invariance 
implies that the boundary measure is the same as the measure on zero sets of holomorphic 
sections used in HU. 

As mentioned above, ‘heat kernel’ random metrics are not like the random metrics of 
Liouville quantum gravity. On a very heuristic level, one may understand the difference by 
thinking of SU(N k )- invariance as a discretization of invariance under the group SDif fu 0 (M) 
of symplectic diffeomorphisms of the background symplectic form (M,ojq). This is far from 
the invariance group of Liouville theory. It is plausible that the only SDiff UJ0 (M)~ invariant 
probability measure on K \^ 0 ] is S Uo . 

This concentration of measure phenomenon, that heat kernel random metrics concentrate 
on the background, is the same phenomenon that occurs for random zeros in [T?i [18] and 
for the quantum Hall point process in [14] . Random zero sets of N zeros or random point 
configurations of N electrons also concentrate at the background metric as N —y 00 . In these 
cases it is customary to dilate the space to obtain configurations of constant density. In the 
case of random metrics, if one dilates small balls of radius ^ around a point z 0 by the factor 

Vk, then the random metrics become metrics on C and the 2-point correlation function of 
the limit measure acquires the new term 5(z — w). 

Acknowledgments. SK is partially supported by the Max Delbriick prize for junior re¬ 
searchers at the University of Cologne, the Humboldt postdoctoral fellowship and the grants 
NSh-1500.2014.2 and RFBR 15-01-04217. SZ is partially supported by NSF grant DMS- 
1206527 and the Bergman trust. We gratefully acknowledge support from the Simons Center 
for Geometry and Physics, Stony Brook University where this paper was completed. 
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2. Bergman metrics 


We briefly review the properties of Bergman metrics, referring to for further back¬ 

ground. As above, (M, u, J) can be any compact Kahler manifold with integral Kahler form. 
The Riemannian metric is g(X, Y ) = a i(JX, Y). The simplest case is that of a Riemann sur¬ 
face, where a Kahler class is the same as a conformal class of metrics with fixed area. Instead 
of parametrizing metrics g = e u go by the Liouville held u we parametrize them by the Kahler 
potential 0, i.e. = cj 0 + idd(p, relative to the reference Kahler form co 0 . 

Bergman metrics of degree k are special Kahler metrics induced by holomorphic embed¬ 
dings 

L s (z) = [s 1 ,...,s Nk }:M^CF N *- 1 

of M into complex projective space. Here, {s^} is a basis of the space H°(M, L k ) of holomor¬ 
phic sections of powers L k of an ample line bundle L —> M with first Chern class Ci(L) = [cuo]. 
Also N k = dim H°(M,L k ). Given a reference basis {s^-} one obtains all others by applying 
an element A G GL(N k , C) to it s A j = ]P AjiSi and induces the embedding 

ts A '■ M —> CF N * 1 , l Sa = A o i s . 

The associated Bergman metric is then, 


I'sa^FS = 


1 N k 

—idd log > 
k 

3 =1 


s j(zj 


(13) 


Since U(N k ) is the isometry group of ujps , the space of metrics is the quotient symmetric 
space V^ k = GL(N k} C)/U(N k ). With no loss of generality one may restrict to SL(N k , C) 
and obtain the quotient SL(N k ,C)/SU(N k ). 

We choose a basis of sections {sj(z)} = (si(z),..., SN k (z)} of H°(M, L k ) which is orthonor¬ 
mal with respect to the reference (background) metric on L k and the corresponding Kahler 
metric ujq = — iiddloghg on M 


1 

V 


n 


Si(z)sj(z)ho ^ = Sij, 


(14) 


where n = dim M. The Bergman kernel of the background metric is the kernel of the 
orthogonal projection onto H°(M, L k ) with respect to the inner product above, and is given 
by 

N k 

B k (z 1 ,z 2 ) = ^Sjizjsjizz) (15) 

3 =1 

Given a positive Hermitian matrix P = Pij the associated Bergman metric is, 


Uab{z) = ^9 a d- b logs i (z)P ij s j (z). 


(16) 


In terms of A G GL(N k , C) above, P = A'A. We introduce the Bergman potential as follows 

(17) 


<t>P = log s i(z) Pij Sj(z) = ^-log\(e A Us(z),Us(z))\ 2 . 


A key property is that /C[ wo ] = Uft. the full space of metrics in a fixed Kahler class 
is the closure of the set of Bergman metrics. Hence /C^j is well approximated by B k for large 
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k, and there are now many results showing that it is well approximated in much stronger 
geometric ways. This approximation problem was raised by S. T. Yau in |2T], see [16] for 
background. 

2.1. Berezin kernel. The key invariant is the Berezin kernel (j5]), given in the above notation 
by 


P 


or in terms of the Bergman kernel 


l (^ i ),^ 2 ))| 2 

K*i)I 2 K*2)| 2 ’ 


p = p k(z 1 , 22 ) 


\B k (z l ,z 2 )\ 2 
B k (zi,zi)B k (z 2 ,z 2 ) ' 


(18) 


(19) 


2.2. Matrix-metric correspondence. The matrix-metric correspondence Eq. (1T61) uses a 
choice of basis {s.,} of H°(M, L k ). Any natural measure on B k must be independent of the 
choice of this basis. We pause to describe such natural measures. 

Any Kahler metric u = oj 0 + idd(p in /C[ Wo ] induces an inner product Hilb fc (</>) on H°(M, L k ) 
by the rule 

(Sl, S 2 )Hilb fc (0) = 

Given a background inner product Go = Hilbfc(0o)> any other inner product has the form 
(si, s 2 )g = ( p gs 1 , s 2 )go where Pg is a positive Hermitian operator on H°(M, L k ) with respect 
to Go- It has a well-defined polar decomposition e A U where U G U{Go) is unitary with 
respect to Go- Its eigenvalues are encoded by the diagonal matrix Aq and its eigenvectors 
are encoded by U. 

In making calculations, we need to parametrize such positive Hermitian operators by pos¬ 
itive Hermitian matrices, which requires a choice of a Go-orthonormal basis of H°(M,L k ). 
Any measure intrinsically defined on the space of positive Hermitian operators will be in¬ 
dependent of the choice of basis. Haar measure and the heat kernel are examples of such 
measures. 


[ si(z)s 2 (z)h k 

Jm n\ 


3. Heat kernel 

In this section we review the heat kernel on 'PaIE Bergman metrics are unchanged if the 
positive Hermitian matrix P is multiplied by a scalar, so we may normalize P so that det P = 
1. Then V N = G/K where G = SL(N, C) and K = SU(N). We denote by t the Lie algebra 
of the maximal compact subgroup K c G and let g = i © p. Let a be a maximal abelian 
subspace of p and let t = dim a. The set of positive roots is denoted by R + . The roots are e, — 
ej, and the positive roots satisfy i < j and have multiplicity one. For SL(N, <C)/SU(N) the 
half sum of the positive roots is the element 8n = ..., ^p-). For background, 

see [I31II2]. 

We refer to the matrix decomposition P = We A U for A = diag(Ai,..., Xn), and U G U(N) 
as ‘polar coordinates’ on V N , where real numbers A j G (— 00 , + 00 ) correspond to the Cartan 
elements of SL(N, C). 

3 In this section we adopt shorthand notation N = N k . 
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The CK (Cartan-Killing) metric is given by 

ds 2 = Tr (P^dP) 2 


( 20 ) 


for P e GL(N,C)/U(N). This metric is bi-invariant under the action of GL(N,C). 

The associated volume form dV on the symmetric space SL(N,C)/SU(N) of positive 
Hermitian matrices with det P = 1 is the bi-invariant Haar measure, 


dV 



N 

A 2 (e A ) d\j 

3 = 1 


m 

[dUu(i)N] 


( 21 ) 


where [dU] is the standard Haar measure on unitary group. 


3.1. Heat kernel measure on SL(N, C)/SU(N). Following Gangolli [10] (Proposition 3.2; 
see also [3], section 2), the heat kernel on SL(N,C)/SU(N) with respect to the standard 
CK (Cartan-Killing) metric is given in ‘polar coordinates’ (A, U) on Vn by 

dfi t = g t (\)dV = C(t,N)^^e-^U^dV. (22) 

A(e A ) 

Here, A(A) = rii<j(-\? — * s th e standard Vandermonde determinant. 

The normalization constant C(t, N) in (1 22]) is fixed by the condition that /z f is the proba¬ 
bility measure f dfi t = 1, 


C(t,N) 


Vn 

27r(\/47rt) 7v2_1 


±N(N 2 - 1 ) 


In deriving this we use the volume of the unitary group 


(23) 


N 

Vol U (N) = (2vr ) N(N+ 1 )/ 2 /Y[j\, 

j= 1 


see e.g. [15]. The factor e i* n ( n2 Pise t H <5jv H 2 and arises because | |<5 tv| | 2 is the bottom of the 
spectrum of the Laplacian. Putting (l2Tj) and (1231) together, we get the following expression 


dut = 


Vn 


27r(\/47rt) Ar2 


-i 


-AJV(JV 2 -1)^ 



N 


A(A)A(e A )e 


4t V 


3 = 1 3 


3 = 1 


[dU] 


[dUjj^y 


(24) 


for the heat kernel measure on SL(N,C)/SU(N) with respect to the CK metric. 


3.2. Geometry of the heat kernel and Brownian motion. In |2J it is proved that the 
mass of the heat kernel concentrates along the exponential image of the U (N)-orbit of the 
5/v-axis in a small annulus centered at 2|( s i A r|f. 

If we write H = diag(A), then the Gaussian factor ^-( Ar “~ 1 )/ 2 i s similar to the heat 
kernel of Euclidean space. But this Gaussian factor must compete with the exponential 
volume growth factor A(e A ) and the factor e~™ N ^ N -P d ue to the existence of a spectral gap 
for A. The well-known factor A(A) pushes the eigenvalues of logP apart. The factor J(H) 

is bounded by e 2 P ,v, P and a simplified expression for the heat kernel is e~ t ^ N ^ +< ' X,pN ^^. 
The maximum of the exponent occurs when A = 2 tSw- 
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Following [2J, let 7 (t) be a positive function with —> 00 as t —> 00 , and let R(t ) be 

a positive function such that R(t)/y/t -A 00 . Consider the annulus 

A(2\5 N \t - R(t), 2\S N \t + R(t)) := {H : 2\5 N \t - R(t) < \H\ < 2\8 N \t + R(t)} C a 

and consider the solid cone 

T (t) = solid cone around the <V axis of angle 7 (i), 

and let 

Q(t) = A(2\5 N \t - R(t ), 2\p\t + R(t)) n T(t). 

Then, according to Theorem 1 of [2], 

I dpt — 7 * 1 , t —^ 00 . (25) 

J U (AT) exp H (t) U (N) 

A Brownian motion proof of this result is given in [3]. It shows that as t —> 00 the mass 
of p t moves off to a component of the ideal boundary (at infinity) of TV. In )|5]we discuss 
this boundary. 

3.3. Scaling and dilation. In the large k limit, the symmetric space metric on Bk, when 
properly scaled, tends to the Mabuchi metric qm on /C[ wo ] (see 0 )- The Mabuchi distance 
function is induced by the Riemannian metric on defined by ||<50||| () = J M (5(p) 2 uj^/n\ 
where = u> 0 + idd(j). We refer to |16j for background. For all k, Bk C /C[ wo ]. If we rescale 

the CK metric gcK,k (ESI) as gi- = e\gcK,ki with = k~ x N k l ^ 2 , then gt- —> Qm on TBk■ Thus, 
a ball of radius one with respect to the usual CK metric gcx,k has radius approximately £*, 
with respect to the Mabuchi distance. It is obviously desirable to consider the heat kernel 
measures for this rescaled sequence of metrics. 

If we rescale the CK metric to gt z = £lgcK,k the corresponding Laplacian scales as A gk —* 
e k 2 A gcKk . It follows that the heat operator scales as 

exp tA gk =expt£ k 2 A gcKk . 

In effect, it is only the time that is rescaled and the rescaled heat kernel is Pk(£ k 2 t, /, P). 

4. One and two point correlation functions of random metrics 

In this section, we calculate the one and two point functions of the random Kahler po¬ 
tential; in the introduction, the latter was stated to be (JUJ) . We use the notation E = 
for the expectation, which at the beginning could be with respect to any V(iVfc)-hivariant 
measure and then specializes to the heat kernel measures. For simplicity of notation we often 
abbreviate N k by N and drop the explicit fc-dependence. 

As mentioned in the introduction, the one and two-point functions are the data required 
to study the mean and variance of the area random variables Xjj. Evidently, 

E k Xu = I E k u, Var(AV) = / E k [uj(z 1 )uj(z 2 )\ - I E^^i^E*, [u{z 2 )\. 

Ju JuxU JuxU 

The integrands are the one- and two-point correlation functions. 
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4.1. One point function of the Kahler potential. The argument here follows PJ. Using 
the integral representation of the logarithm 

, 1 

log a = - + 7 — 
r 

where 7 is the Euler constant, we can rewrite the expectation value of the Kahler potential 
<t> P m relative to the Kahler potential 0/ of the background Bergman metric as 

E k[4>p(z) ~ <Mu)] = E k 

— y lim — + 7 — f x T ~ x dx [ e- TveAuwt dfi u (27) 

r J o J Vn 

where introduced the matrix T ji = xsj{z)si(z)/\s{z)\ 2 . The integration over the unitary 
group can be carried out using Harish-Chandra-Itzykson-Zuber formula, see e. g. [T51 122] for 
background. Namely, for any two Hermitian matrices A and B with eigenvalues aj and bj 


1 , s(z)U'e A U s(z) 

-log- 


k 




X 


T — 1 „ — OLX 


e~ ax dx + 0 (t), 


(26) 



[dU] 

Vol U(N) 


exp (fiTrAUBrf) 



N(N- 11/2 ^ 

A (a) A (6) 


(28) 


It is not hard to check that this expression is well defined even if some of the eigenvalues 
coincide. This is the case for the matrix T which has N — 1 zero eigenvalues and one non-zero 
eigenvalue equal to x. Hence the integral on the right hand side of (|27j) is z -independent 
and we immediately conclude that the expectation value of the Bergman metric is equal to 
the background Bergman metric 

In fact, this is true for any eigenvalue-type measure P, since the HCIZ integral depends on 
eigenvalues only. Note, that so far no assumptions on k have been made in this calculation. 
Considering now the limit of k large, we invoke the Bergman kernel expansion [9] to show 
that the background Bergman metric tends to the reference Kahler metric, 

kfyj, 06 = U 0ab + 0 (l/k). 

Now we would like to consider the variance of Xy from its mean ojq(U). 


4.2. The two point function. In this section we prove the formula (PJ) for the two point 
function, and then discuss its asymptotics in various regimes. In the terminology of Cl we 
are finding a ‘bi-potential’ for the variance. Although the calculations of this bi-potential are 
completely different from the case of random holomorphic sections in d, the final formulae 
are somewhat similar and when t 00 they are identical. 

Recall that the two-point function (JJJ) of the Kahler potential is the sum of the background 
term cj) I (z)(pi(w) plus the variance term 7 2 ,fc(p)- Instead of using the formula (( 8 ]) for this term, 
we take the approach of writing 


K 2 , k (z,w) - 0/(z)0 7 (?/;) = — E fc 


log sjz^We^sjzj) lQg l{z 2 )We A Us(z 2 ) 


K*)l s 


1 ^ 2 ) p 


= — lim (l 2 k(t, P,Ti,t 2 ) + p-independent terms) 

k Z ri,r 2 ^0 v ’ ' 


( 29 ) 
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where we do not write down the p-independent terms, since ultimately we are interested in 
the dependence of the correlation function on coordinates, which enter only through p(z\, z 2 ). 
Here 

p poo p 

P) H, t 2 ) = // x T 1 1 ~ 1 xl 2 ~ 1 dx 1 dx 2 e~ TreAuiS,ut dp t . (30) 

JJo Jv N 

and we introduced the matrix + x 2 3 ~-tj p 2 ^ ^ l ias ran k 2 with two 

non-zero eigenvalues given by 

0i,2 = i (xi + x 2 ± V (1 - P){%i - x 2 ) 2 + p(x 1 + x 2 ) 2 j . 

Applying the HCIZ formula (|28|) to the unitary integration in Eq. (130|) we obtain 


D /2 N!(N-iy. _ e-^ 

y ’ (**) N - 2 (A-^)nf, 2 (e il -e A J)n;U(e A “-e i <) 

where we used the fact that the integration measure is symmetric in eigenvalues A’s. 

Now we use the explicit form of the heat kernel measures p t (t24|) . The integral I 2) k(t, p, Ti, t 2 ) 
in Eq. (129|) with this eigenvalue measure can be written as 


/ 


k(t,p,Ti,T 2 ) = {- 1 ) N{N 1 )/2 C'(t, N) 

N ( 

\d\j A(A)Ai 2 (e A ) exp I - 

3 =1 V 



Vo\U(N) f[°° x^xT'dx.dx, 

(2 n) N [ J \/./ 0 (0 1 0 2 )^ v - 2 (0 1 -0 2 ) 

N N \ 

a j + iy x i ~ ^ ieAl _ <^ eA2 ) > 

3 =1 3 =1 / 


where we defined the partial Vandermonde determinant Ai 2 (e A ) = ri 3 <j<Kiv( eAj ~ eAi )> which 
excludes the first two eigenvalues e Al and e Az . The ^/-integration enforces the delta-function 
constraint in the measure (1241) . 

Using antisymmetry of A (A) under exchange of two eigenvalues, A 12 (e A ) can be replaced 
by (—l) 1+A A Ar “ 1 )/ 2 (jV — 2)! e^h=3^ -3 ) Ai inside the integral, which leads to the further simpli¬ 
fication 


h,k(t, P, A, r 2 ) 


Y ° W ^ C{t, N)N\(N - 1)\{N - 2)! 


f C°° x ^ 1 l x T 2 1 dx\dx 2 
JJo ( 0102) 7V_2 (01 - 02 ) 





(31) 


Thus after the HCIZ integration we got rid of most difficult factor A(e A ) and left with a 
Gaussian integral with a polynomial amplitude, except for the terms 0ie Al + 0 2 e Aa in the 
exponent. Note that the P-dependence is entirely in the variable p inside 0i, 0 2 - 
The next step is to calculate the integrals in the second line of f|3ip . 


01 , 0 2 ) 





(32) 
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Our strategy is to Taylor-expand the exponent of (152]) in powers of e A , 


0 -<t> ie A i—</> 2 e A 2 _ 


E 

m,n =0 


_ nh+mX , 

, , & 1 


n\m\ 


then perform gaussian integration in A term-by-term, and finally re-sum the resulting series, 
i. e. un-do the Taylor expansion. We use the identity 


N 

n 

j = 1 

to compute the eigenvalue integral and get 

(n — m) 


d\j ) A(A) = (2ir) w/2 (2i) wl/2 A( fl ) e'V'-M 


(33) 


N -3 


01 , fa) — 


n,m =0 
poo 


n\m\ 


-*)"(-«”* Il("- 0(">-0 


1=0 

2 1 4 v^iV-3|f, , -,n2 


dy e t(.n+iy} 2 +t{m+iy) 2 +tJ2i =0 (i+*») _ 


(34) 


Due to the factor riil: 0 3 ( n — 0( m — 0> terms with m,n < N — 2 have coefficient zero, so 
we can shift summation indices n —» n — (iV — 2), m —> m — (iV — 2). Integrating over y in 
and plugging the result back to (15T|) we obtain 


r u „ ^ ^ \ t(lV— i) 2 /iv f f x i x 2 dx\dx 2 

hAt,P,n, r 2 ) = -e yy o — 

n!m! 


n,m=0 


Now we can re-sum the series using the identity 


1 

2nt 


dAidA 2 e 2t 


-A.(A 2 +A 2 )-^ie aAl+A 2-0 2 e aAl - A 2 _ 


E 

n,m =0 


-01 ) W ( 02) m e | a 2( n+m )2 + t ( n _ m) 2 


n\m\ 




Replacing the series in (]35l) by the integral, and changing variables X\ —> e “"' i 0 n x\ 1 
X 2 —> e _aAl_f_ n X 2 makes it possible to carry out the gaussian integration in Ai, giving 


02 ,k(t, p, Ti, T 2 ) — 


o~t /2 


\j2ixt 


00 x\ l 1 xl 2 2 l dx\dx 2 


/•LXJ 

(0n9 01 - 0 2 c0 2 ) / dA 


Interestingly, the factor coming from the spectral gap has now disappeared from 

the formula. It follows from the integral representation Eq. (I26]h that the singular in 7T >2 
terms in (129|) are p-independent. Therefore after taking the derivative of / 2 ,fc(U p, t±, r 2 ) with 
respect to p, we can set T\ = r 2 = 0. Since the measures p t depend on t we henceforth denote 
the corresponding term in the two point function by J 2j fc(f, p). Using 


d p (j) 1,2 — 


xix 2 


d, 


2xix 2 


01 02 ^ 01 -02 (01 - 02) 3 
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we get 


d p h,k{t,p) 


2e~ t/2 

\Z2nt 
sinh A 


dxidx 2 


dXe-^ 2 -^-^- 


■( 0 i + 02 + ( 0 i - 0 2 )( 0 ie A - 0 2 e A )) 


(0i - 02) 3 

Now we introduce new coordinates (r, 6) as 

r cos 6 = s/p(x\ + x 2 ), r sin 6 = y/l — p (x 2 — xi), 

with the range 

r G [0,oo), 6 1 G [—cp +a] where cos a = \fp. 

In terms of r and 9 the integrals can be written as 


d P h,k{t,p) = 


=,-*/2 


V2nt sjp{l - p) Jo 
cos 6 1 


rdr 


dO 


cos 


o + Vp 


cos 


e — 


dX e 


0-VP 


— A- A 2 —r ( 22^ cosh A+sinh A j 


sinh A. 


_r 2 yp 2r V yp 
Introducing new variable Xy/o = cos 0 and rearranging exponents, we get 


d p I 2>k (t,p) = 


o-t/2 


y/2irt -\/l — p 


rdr 


-i /Vp 


dx 


\/l - px 2 J-c 


dX e 2 1 


-Aa 2 +a 


—rx cosh A 


£(, 

V 


—r sinh A 


_r sinh A 


sinh A 


) +-((z + 1> 


—r sinh A 


+ (x- l)e rsinhA ) 


Integrating over r we obtain 


d p h,k{t,p) = 


-m 


■i/Vp 


dx 


a/27 Tt a/1 — p J 1 
x log 


n/T" 


<ZAe 2U 


/0X Z J —oo 


x cosh A — sinh A A „ . . . x 2 cosh A — sinh A 
j + 2 sinh A 


_ x cosh A + sinh A J x 2 cosh 2 A — sinh 2 A _ 

After integrating the log term by parts, we can perform the x integration 

n/v/p (] x e~Jt x2 cosh A sinh 2 A 


Q r „ .a 2 1 2e~ t ' 2 a/W 

dpl 2> k[t, p) — r7T~: 

p \J2irt p 


dX 


21 
P 


o-t /2 


V 1 ~ P 


a/ 27rt p 


dXe 2 t' 


cosh A 


log 


px 


2 x 2 cosh 2 A — sinh 2 A 


y/ coth 2 A — p + a/1 — p 


(36) 


a/ coth 2 A — p y/ coth 2 A — p — a/1 — p 

This completes the calculation of (J9]). 

An important application of these asymptotics is to calculate the variances of the linear 
statistics X v (J3J) and its smooth analogue Xf. This can be done precisely as in Section 4 of 
m and Section 3 of [TS], but substituting the formula for I 2>k (t,p) for Q k . The details are 
rather lengthy and will be presented elsewhere. 

The remainder of the article is devoted to the asymptotics of I 2)k {t, p) as k — » oo with 
various relations between t and k. To understand the various regimes, it should be recalled 
that the metric rescaling in H3.3l is necessary to ensure that balls in the Cartan-Killing metric 
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on Bk maintain their size with respect to the limiting Mabuchi metric on /Cp, 0 ]. Without 
rescaling, the Brownian motion relative to the CK metric is probing metrics only at distance 
dk — £&£ from the initial background metric with respect to the limiting metric. 

In all of the regimes, the key to finding the scaling asymptotics is to work out the behavior 
of 


A(t : p) : = 


log 


a/CO tll 2 t — p + y/1 — p 


(37) 


\J cotli 2 1 — p \/cotli 2 1 — p — a /1 — p 

as k —y oo, where t may depend on k. As will be seen below, the factors of (a/27t t)~ 1 e~ t / 2 
in front of the integral are always cancelled, leaving the prefactor ■ y -S. By (1431) . for pairs 


(z,z + -^=), we have p —* e l“l 2 , and Ak(t,p) has a limit, which depends on whether or not 
we also send t —» oo. 


4.3. The limit as k —> oo for fixed t. This regime corresponds to letting the Brownian 
motion with respect to the Cartan-Killing metric evolve for a time t, and as discussed in §3.31 
the ball of radius t in fjcK.k metric is shrinking in size with respect to the Mabuchi metric 
and has c4-radius equal to ept. 

The asymptotics in this regime could in principle be derived from the formula ()8j) of [9]. 
But more explicitly, we note that p —* 0 as k —> oo off the diagonal. Expanding at small p, 
we get 


a/27 rt 


,- t /2 


dXe 


cosh A ,„.A/coth + 


and the first term here cancels the first term in 
(zi, Z 2 ) fixed then p —> 0, we get 


a /coth 2 A — p \/ cotli 2 A — p — a/1 — p 

Thus in the regime when we 


(38) 

hold 


h,k{t,p) — ao (t) + a\{t)p H-, 


where ao, ai,... are constants independent of p. To obtain the 2-point correlation function of 
the Kahler metric, we then take four more derivatives. The constant a 0 does not contribute 
to the answer and we see that the two point correlation function is the free background term 
kfyi (^i)^/( 22 ) plus a term exponentially decaying off the diagonal like C 2 (t)e~ kDl ^ Zl,Z2 \ 


4.4. The limit as t —> 00 for fixed k. Now we apply steepest descent to the second integral 
as t —> 00 , and keeping k, Nk fixed. We obtain 

<V 2 ,fc(oo,p) := lim d p I 2 , k (t,p ) = lim — - ~(2t + log(l - p) + 0(l/t)) = -— • 

t-*oo t->oO pp p 

Thus we have, 

h,k{oo,p) = Li 2 (p). (39) 

As mentioned above, this is the same formula as (1T2]) . In the scaling limit around the 
diagonal with pairs of points of the form (z, z + A=) we have p ~ e~\ u \ . In the next section 
we connect this limit with the correlations between zeros of Gaussian random holomorphic 
sections, see (1551) for a more precise statement. 
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4.5. The metric scaling limit with t —$■ te k 2 . The goal now is to evaluate I 2 
asympotically as k —)■ 00 . This scaling keeps the (4-balls of uniform size as k — > 00 with 
respect to the limit Mabuchi metric. Thus, as k changes the Brownian motion with respect 
to (jf ; probes distances of size t from the initial metric ta 0 for all k. 

In order to apply steepest descent for A G [0, 00 ] we change variables A —> £ k 2 \ in so 
that the exponent becomes 



The saddle point occurs at A — t, and the critical value of the exponent is |. The prefactor 
e -b 2 j n front of the integral, with t —> £ 2 k t again cancels the critical value of the phase and 
the singular term. The integrand in (136|) at the saddle point is asymptotic to 


coth 2 (e k 2 t) — p J coth 2 (e k 2 t ) - p - y/1 - p 


log 


\]coth 2 (e k 2 t) -p + y/l-p 


2£ k 2 t 
a _|_ 


2 / _—2^ 


V 1 - p V 1 - p 


log(l - p) + 


(40) 


so 


exactly as in (l39]h 


d p l2,k{e k 2 t )P ) ~ 


log(l - p) 

? 

P 


4.6. Spatial scaling. In this section, we consider the scaling asymptotics of / 2 ,fc(p) discussed 
in §1.11 The only new element in the calculation is the scaling asymptotics of the Berezin 
kernel p. The calculation of h,k(p) above does not change. 

The main input into the scaling asymptotics is the following facts about the Berezin kernel 
(see nzmsi for background and references). Off the diagonal one has 


p ~ e 


-kDj(zi,Z2) 


(41) 


where the Calabi diastasis function is given by 


Di(zi, z 2 ) = 24 ) + </>i(z 2 , z 2 ) - 0/(^1, z 2 ) - 4 >i(z 2 , zi). (42) 

Here 4>i(z,z) is the Kahler potential (|I7|) of the background Bergman metric ccorre¬ 
sponding to the identity matrix P = /, and (pi(zi,z 2 ) is its off-diagonal analytic extension. 
The Berezin kernel has a scaling limit on the ‘near-diagonal’ where the distance d(zi,z 2 ) 

between 24 and z 2 satisfies an upper bound d(zi, z 2 ) < b ( b G M + ). By comparison, 

1 /2 

p —> 0 rapidly on the ‘far off-diagonal’ where d(z\, z 2 ) > b (^|^) , in the sense that for all 

6, R > 0, 

W j Pk(zi, z 2 ) = 0(k~ R ) uniformly for d(zi,z 2 ) > b 

Here, V' ? stands for the j-th covariant derivative. 

The scaling asymptotics of this kernel near the diagonal may be described as follows. Let 
2 G M. Then 



Pk 


z + 7I’ z + 7t) =e "A 1 + R *( u ’ v )1 • 


(43) 
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where 

\R k (u,v)\ < % \u-v\ 2 k-R 2+E , |Vi2 fc (u)| < C 2 \u-v\k~R 2+£ , 

\ViR k (u,v)\ KCjk- 1 / 2 ^ j> 2 , 

for |w| + |u| < fr-y/log k. 

It follows that the scaling asymptotics of the variance term for random metrics is given by 

h,k{£l 2 t,P k (z,z + u/Vk)) ~ Li 2 (e _|u|2 ), (44) 

just as in the limit as t —> oo first, considered in H4.41 

4.7. Energy entropy scaling. Another natural scaling comes directly from the density 
(124)1 . We separate out the ‘action’ from the ‘amplitude’ and express it in terms of the 
empirical measures d/i\ = A- Y^j =i of the eigenvalues of A 

A(A)A(e A )e _ 4t ^' =1 = 

= e N k(f M /M lo sh“d^AbO^A( 2 /)+(/ M f M i°g\e x -e y \ d v\(x) d fj'k{y)))+^r f M x2d ^ . ( 45 ) 

To give all terms the same order in N k we need to rescale the time to t —> t/N k . Then the 
scaled measures p k satisfy a large deviations principle with the rate function given by the 
exponent of (l45j) . 

5. The geodesic boundary of B k and configurations of zeros 

In this section, we explain how the result of 1 14.41 is essentially the same as the theory 
of zeros of random holomorphic sections. As t —>• 00 , the mass of the heat kernel gets 
concentrated on a part of the boimdary of V^ k corresponding to ‘singular metrics’ given by 
zero sets of holomorphic sections. 

Symmetric spaces of non-compact type have several different notions of boundary and 
several types of compactihcations. For background we refer to HU The boundary relevant 
to the heat kernel measures and their t —> 00 limit is best stated in terms of the Bergman 
metrics themselves and their limits along geodesic rays of V^ k . 

Definition 1. The weak* compactification of B k is B k U dB k where dB k is the set of limit 
points (i.e. endpoints) of the Bergman metrics along Bergman geodesic rays co k {s). 

In fact, the only relevant boundary points are the ones arising from the geodesic ray 
starting at the balanced metric x k ~ 0 and with initial velocity in the direction of S jv, 
together with the endpoints of the U(N k )- orbit of this ray. 

A geodesic ray in the space of Bergman potentials is a one-parameter family of metrics 
whose potentials have the form, 

A=(: lo «E e “ , l s j'WI 2 - 

3 

where in SL(N k ,C)/SU(N k ) the ray starts at the origin and has initial vector (U,A). We 
note that 

^\ogJ2e tX ^s^(z)\ 2 = t -^ + ^\og(\s^ ax (z)\ 2 + Y, e t(A ^ Amax) I^W| 2 V (46) 

j \ max / 
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Here, A max is the largest of the A j and s^x is the corresponding section. Clearly, 

\ log Ej etXj 0) I 2 - sup M l log E j etXj I sf (z) | 2 

= \ log (kExWI 2 + EA^max e * (Ai_W) | S 7(^)| 2 ) (47) 


-SU Pm ±log |sEx(-)| 2 + EA^max 6 '^ Amax) | sf (z) | 2 ) . 

This family of potentials is bounded above by 0 and is pre-compact in L P (M ) for all 1 < 
p < oo. The sup in the second term is itself bounded above by and therefore we may 
remove the sup without changing the limit (i.e. the sup was only needed to get rid of the 
tA ™ ax term). We then observe that for any Bergman geodesic ray, 


lim 

t—} OO 



JJ 


E 

Aj^max 


0 t(Xj Amax) | gU ^ 



(*) 


—> 0. (48) 

V-(M) 


Indeed, 

Ft(z) ■= ( !<&*(*) I" + J] e *(Aj-A m ax)| s U( z )| 2 J (49) 

is monotonically decreasing to (sEx^)! 2 f° r each z. Therefore its logarithm monotonically 
decreases to log |sE x (^)| 2 - If we subtract | log |s max (^)| 2 then the difference is always > 0 
and we may remove the absolute values, and then apply the monotone convergence theorem 
to take the limit under the integration sign. But the limit equals zero almost everywhere 
since each term tends to zero, proving ()48l) . 

We now consider the (1, 0) forms obtained by taking d of the potentials. Taking d kills 
the * A ™ ax term, and gives, 


Eiog£ e ‘y s p*)i 2 = imo g i4L*wr+ E e ,(A '- x -“>| s ywr . (so) 

j y Aj^max J 

Proposition 5.1. The weak limits of Bergman metrics along geodesic rays are generically 
given by the normalized zero distributions of holomorphic sections of L k , i.e. as t —>• oo, 

rddlogY e tAj |sf (z)\ 2 —> i(?91og|s| 2 . 
j 

If the highest weight has multiplicity r, then the limit is |<9<91ogEj=i I s / 2 where 1 < r < n 
and {sj}t = i is any set of sections in H°(M, L k ). 

Proof. We take dd of the two terms of ()48l) to obtain, 


-<9<91og Ee fAj I^WP 
k J 


-A —dd log 
k 


(*) 


2 


(51) 
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in the sense of distributions. That is, if we integrate against a smooth (n — 1, n — 1) form 0 
(i.e. a smooth function if the dimension of M equals one), then 

/ M 0 A \dd log I 2 + E^max e^-^l sf (z) I 2 : 

= I M 99 ^ A l log (|sEx(-)| L> + EA^max et(Aj_Amax) kj / (^)| 2 ) 

and by (|48l) the right side tends f M ddi /> A | log |s ^ av (z)| 2 . 

Multiplicity r means that A max has multiplicity r, and then one sums over the associated 
sections Sj. □ 

A key point for this article is that each of the weights in 5^ is distinct. Hence r = 1 in that 
direction, and therefore the boundary points in the direction Sn and its U( Retranslates all 
consist of 5 -functions tBB log I s(z)| 2 of holomorphic sections of L k . Every section arises in 
the U (iVfc)-orbit. 

5.1. Relation to zeros of holomorphic sections. In mi, the authors found a bi¬ 
potential Qn € L l {M x M) for the pair correlation function of zeros of Gaussian random 
holomorphic sections 

1 n 

= ^e' |c|2 dc, s = c o s i > ( 52 ) 

i=i 

on L k , where {s^} is an orthonormal basis and dc is 2n- dimensional Lebesgue measure. This 
Gaussian is characterized by the property that the 2 n real variables Re Cj , Im c.j (j — 1., n) 
are independent Gaussian random variables with mean 0 and variance A i.e., 

Ecj = 0, E CjCk = 0, E CjCk = Sjk ■ 

The current of integration Z s over the zeros of one section is given by the Poincare-Lelong 
formula 

= ^^ddlog\f\ = ^^^loglls^H +0/(2r) , (53) 

7T 7r 

As mentioned above, the limit as t —> oo of random Bergman metrics along the rays above 
must be random singular metrics, and we claim that the limit ensemble is equivalent to the 
Gaussian one. Indeed, the limit measure is U(Nk)- invariant and there exists just one such 
measure up to equivalence, namely the Gaussian measure above. This explains why the 
bipotential 47 t 2 Qk of is the same as the t —> oo limit of the two-point, function l 2 ,k{p) 
of the heat kernel ensemble in 1 14.41 corrobrating that this limit ensemble is that of random 
zeros of sections. 

5.2. Smooth linear statistics and asymptotic normality. Similar to the area random 
variable Xu but easier to work with, is the smooth linear statistic 
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where / G C°°(M ) is a smooth test function. It has a much smaller variance than Xu since 
the effect of ‘zeros along the boundary’ is smoothed out. In [18] it is proved that 


Var (X f ) = k~ l 


'03) 

16(77 


IIA/II*+ 0(*rK«) 


(55) 


It was proved by Sodin-Tsirelson [19] in certain cases and then m for the present setting 
that the random variables Xf obey an asymptotic central limit theorem. Namely, if we 
normalize Xf to have mean zero and variance one, then 


X f - E (X f ) 
VVar( X,) 


(56) 


tends in distribution to the standard Gaussian distribution Af(0, 1) as k —>• oo. That is, 


kl/2 ( x ' ~ l L foJ ") ^M 0 . v^riiwib) 

in the weak sense of convergence of distributions as k —> oo, where is a certain positive 
universal constant and A/"(0, a) denotes the (real) Gaussian distribution of mean zero and 
variance a 2 . 

The same results hold in the limiting case of the heat kernel measures when t —» oo. It 
would be interesting to investigate the analogous variance and asymptotic normality results 
for heat kernel measures in the other regimes as k —> oo, t —> oo. The formulae of this article 
combined with the techniques of [1?, [18] make this possible. 
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